% 中断概率与SNR或R的关系图

% 生成数据
SNR_dB = -10:5:30;

for I_SNR = 1: length(SNR_dB)
    rho = 10^(SNR_dB(I_SNR)/10);  
    r=1;
    Out_SNR_1(I_SNR) = 1 - exp((1-2^r)/rho); % r=1中断概率计算公式
    r=3;
    Out_SNR_3(I_SNR) = 1 - exp((1-2^r)/rho); % r=3中断概率计算公式
    r=5;
    Out_SNR_5(I_SNR) = 1 - exp((1-2^r)/rho); % r=5中断概率计算公式

end

% 绘制图形
subplot(2, 1, 1); % 第一个子图
hold on
plot(SNR_dB, Out_SNR_1, 'ro-');
plot(SNR_dB, Out_SNR_3, 'bo-');
plot(SNR_dB, Out_SNR_5, 'co-');
legend('Out(r=1)','Out(r=3)','Out(r=5)', 'Location', 'northeast'),

xlabel('SNR');
ylabel('中断概率');
title('中断概率与SNR的关系');
% 添加网格线
grid on;

% 绘制第二个子图

% 生成数据
R = -1:1:10;

rho5 = 10^(5/10);
rho10 = 10^(10/10); 
rho20 = 10^(20/10); 

for Number = 1: length(R)
    
    %信噪比为5db,当db值转化为线性数值
    Out_5(Number) = 1 - exp((1-2^R(Number))/rho5); % SNR=5中断概率计算公式
    
    %信噪比为10db,当db值转化为线性数值
    Out_10(Number) = 1 - exp((1-2^R(Number))/rho10); % SNR=10中断概率计算公式

    %信噪比为20db,当db值转化为线性数值
    Out_20(Number) = 1 - exp((1-2^R(Number))/rho20); % SNR=20中断概率计算公式

end


% 绘制图形
subplot(2, 1, 2);  % 第二个子图
hold on
plot(R, Out_5, 'ro-');
plot(R, Out_10, 'bo-');
plot(R, Out_20, 'co-');
legend('Out(SNR=5db)','Out(SNR=10db)','Out(SNR=20db)', 'Location', 'northwest'),

xlabel('传输速率 R');
ylabel('中断概率');
title('中断概率与传输速率R的关系');

% 添加网格线
grid on;
